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Abstract 

Plant growth and productivity are greatly affected by drought, which is likely to become more threatening with the 
predicted global temperature increase. Understanding the genetic architecture of complex quantitative traits and 
their interaction with water availability may lead to improved crop adaptation to a wide range of environments. Here, 
the genetic basis of 20 physiological and morphological traits is explored by describing plant performance and growth 
in a Brassica rapa recombinant inbred line (RIL) population grown on a sandy substrate supplemented with nutrient 
solution, under control and drought conditions. Altogether, 54 quantitative trait loci (QTL) were identified, of which 
many colocated in 11 QTL clusters. Seventeen QTL showed significant QTL-environment interaction (QxE), indicating 
genetic variation for phenotypic plasticity. Of the measured traits, only hypocotyl length did not show significant gen- 
otype-environment interaction (GxE) in both environments in all experiments. Correlation analysis showed that, in the 
control environment, stomatal conductance was positively correlated with total leaf dry weight (DW) and aboveground 
DW, whereas in the drought environment, stomatal conductance showed a significant negative correlation with total 
leaf DW and aboveground DW. This correlation was explained by antagonistic fitness effects in the drought environ- 
ment, controlled by a QTL cluster on chromosome A7. These results demonstrate that QxE is an important compo- 
nent of the genetic variance and can play a great role in improving drought tolerance in future breeding programmes. 

Keywords: Antagonistic fitness effect, Brassica rapa, drought, genotype-environment interaction, plasticity, root/shoot ratio, 
stomatal conductance. 



Introduction 

Plant growth is greatly affected by environmental abiotic 
stresses, of which drought is the most common factor imped- 
ing crop productivity. Drought is likely to become more 
threatening with the predicted global temperature increase 
(Smith and De Smet, 2012). Three categories of plant adap- 
tive strategies to drought have been recognized: drought 



escape by early flowering, drought tolerance via increasing 
water use efficiency and drought avoidance via reduced tran- 
spiration and increasing water uptake (Levitt, 1972). 

Evaluating those responses in many genotypes in sev- 
eral environments may show phenotypic plasticity, which is 
defined as the ability of an individual organism to alter its 
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physiology/morphology in response to changes in environ- 
mental conditions (Schlichting, 1986). When this plasticity 
differs between genotypes (i.e. when there is genetic variation 
for it), it is classified as genotype-environment interaction 
(GxE) (Via and Lande, 1985). Understanding GxEbetterwill 
provide a solid foundation for genetic improvement of stable 
crop productivity and will help to identify superior and stable 
alleles/genotypes across different environments (Zhang et al. , 
2010). The genetic basis of the observed GxE can be identi- 
fied by genetically dissecting plant physiological and morpho- 
logical responses to environments via quantitative trait loci 
(QTL). This specifies the genetic component of GxE and is 
expressed as QTL-environment interaction (Q X E) (Malosetti 
et al, 2004; Boer et al, 2007; Tardieu, 2013). Different QTL 
effects can occur if the allele underlying the QTL is strongly 
expressed in one environment but weakly in another, or if the 
allele has opposite effects on the same trait in different envi- 
ronments (Mackay, 2001; Sukhwinder et al, 2012). A QTL 
for which one allele has opposite (pleiotropic) effects on the 
phenotype in two different environments can lead to fitness 
trade offs, elevating fitness in one environment but depress- 
ing it in the other environment. Trade offs can be maintained 
in nature (e.g. by antagonistic pleiotropy), when alleles at a 
locus underlying a fitness component show clear home-site 
advantages (Rose, 1982; Anderson et al, 2013). Therefore, 
considering such antagonistic fitness effects is crucial while 
selecting for desirable QTL during marker-assisted breeding 
programmes. 

To facilitate improving marker-assisted breeding pro- 
grammes, a model crop plant is required. The Brassica genus 
has the smallest genome size, the complete genome sequence 
of Brassica rapa (Wang et al, 2011), close relationships with 
the plant model species Arabidopsis thaliana and genome 
analysis tools, provided in the Brassica database (BRAD) 
(Cheng et al, 2011), so B. rapa is a useful dicot model crop 
for genetic and comparative studies. 

The present study focused on drought avoidance, which 
enables plants to maintain a high fitness level in drought con- 
ditions. Therefore, QxE on growth-related traits were investi- 
gated in a B. rapa recombinant inbred line (RIL) population 
grown on a sandy substrate under control and drought envi- 
ronments. This work identified several QTL for main effects 
and QxE and found an antagonistic fitness effect for a sto- 
matal conductance/shoot biomass QTL, with the same allele 
reducing stomatal conductance under drought and increas- 
ing it under normal watering conditions, while contributing 
to higher shoot biomass in both environments. 



Materials and methods 

Plant material and experimental setup 

The RIL population (F7) used here was previously developed by this 
study group from a cross between a Yellow Sarson (R-o-18) (d 1 ) and 
a Caixin type (L58) (9) and genotyped with 270 markers (Bagheri 
et al., 2012). The RIL population was screened three times under 
control (continuous watering for 3 weeks) and drought (normal 
watering for 1 week, then plants were left to dry out) environments. 
In all screens, plants were grown in 13-cm-deep square black plastic 



pots. Each pot was filled with 1.5 kg dried river sand and all pots 
were watered until saturation with 1100 ml nutrient solution (1, 1.1, 
5.9 mmol T 1 , N, P, and K respectively). The same nutrient solution 
was used for watering plants every 2 days. Two seeds were sown per 
pot and 4 days after germination, seedlings were thinned to one per 
pot. Seven days after germination, watering was withheld as drought 
treatment, while the control treatment was continuously watered. 

Initially, a pilot experiment was performed using 30 randomly 
selected RILs and both parental lines, with three replications per 
genotype per environment, to test if the drought treatment would 
reveal significant differences between RILs and between the two envi- 
ronments regarding total leaf fresh and dry weight. Subsequently, a 
full RIL screening experiment was performed in which 140 RILs 
and both parents were phenotyped for the 20 studied traits under 
both environments with three replications per RIL and six replica- 
tions per parental line per treatment. Finally a QTL reproducibil- 
ity experiment was performed to confirm the different phenotypes 
for contrasting alleles at four identified QTL by screening 27 RILs 
selected for their discriminating genotypes, with three replicates per 
RIL per environment. 

All experiments were carried out under controlled greenhouse 
conditions under al6/8 light/dark cycle (22.3/20.3°C, mean relative 
humidity 77.8/81.3%). The experimental setup involved a complete 
randomized block design with one plant per RIL and two replicates 
for each parent per block. 

Plant phenotyping 

In the full RIL screening and QTL reproducibility experiments, 20 
traits were analysed under control and drought environments. These 
traits were chosen as the ones describing as best as possible the dif- 
ferent aspects of plant performance. Directly before harvesting, 
when less than 5% of plants had visible flower primordia, the num- 
ber of leaves was counted. Chlorophyll content was measured (only 
in the full RIL screening experiment) using a SPAD-502 chlorophyll 
meter (Minolta, Japan). For this measurement the average of three 
leaves per plant per replication per treatment was taken. Leaf sto- 
matal conductance was measured using a leaf porometer (Decagon 
Devices, USA) for one fully expanded leaf per plant per replication 
(either the 3rd or 4th leaf). Thereafter, total leaf fresh weight (LFW) 
and dry weight (LDW) and the dry weight of the 3rd and 4th (i.e. 
fully expanded) leaves (3,4DW) was measured. Dry weights were 
determined after drying plant materials at 65°C for 4—5 days until 
weight constancy. 

Leaf area (LA) of the 3rd and 4th leaves was measured using a 
Licor LI-3100 (Licor, Lincoln, NE, USA), and subsequently their 
combined specific leaf area (SLA) was calculated as LA divided by 
3,4DW, as well as the dry weight ratio between 3,4DW and LDW. 
Hypocotyl length was measured using a ruler, and hypocotyl DW 
(HDW) was determined. The shoot DW (SDW) was calculated as 
the sum of LDW and HDW. 

Subsequently, root systems were washed carefully to remove 
adhering sand, placed in a plastic tray filled with water, spread 
and scanned with a flatbed scanner. From this, the total root sys- 
tem length (RL), root volume (RV), and root diameter (RD) were 
measured using WinRhizo (Regent Instruments, Quebec, Canada). 
This was used to calculate the RL-to-SDW ratio (RL/SDW), which 
illustrates the aboveground matter that is supported by a given RL. 

Thereafter, roots were dried to measure root DW (RDW) and to 
calculate the root-to-shoot DW ratio (R/S). Similarly, to indicate the 
relative investment in shoots or roots, the shoot-to-total plant (shoot 
+ root) DW ratio was calculated (S/SR), for which total plant DW 
was calculated as the sum of SDW and RDW. Finally, the leaf water 
content (LWC) was calculated as (LFW - LDW) / LDW. 

Statistical and quantitative trait loci analysis 

Statistical analysis was performed on raw data of each experiment 
using GenStat for Windows 1 5th edition ( VSN International, Hemel 
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Hempstead, UK). Analysis of variance (ANOVA) was used to test 
the significance difference between treatments, lines, and interaction 
(GxE). Heritability was estimated as implemented in GenStat. In 
the linear mixed model, genotypes were fitted as random and blocks 
as fixed. The generalized heritability measure used, as described by 
Cullis et al. (2006), and in a more general context by Welham et al. 
(2010), is given by: 

mean(pev(g ; )) 
h = 1 x — - 

where the set of predicted genotype means (Best Linear Unbiased 
Predictors) are g l ... g N with prediction error variance pevfc) and 
estimated genetics variance component cr g 2 . Pearson correlations 
were calculated using GenStat. 

Data from the 20 traits analysed in the full RIL screening experi- 
ment were used for QTL mapping using a multienvironment analysis 
(ME A) approach, which accounts for GxE, as implemented in the 
QTL library in GenStat. A step size of 10 cM, a minimum cofactor 
proximity of 50 cM, a minimum separation of selected QTL of 30 
cM, and a threshold of -loglOP = 2.8 were used for QTL analysis. 
Following the mixed-model approach described by (Malosetti et al. , 
2004; Boer et al. , 2007), first the whole genome was scanned using 
simple interval mapping and then, based on that, cofactors were 
selected for two rounds of composite interval mapping. Thereafter, 
a final QTL model was selected using backward selection on the 
selected cofactors, where it estimated the allelic effect of each of 
QTL in each environment, the effect of QxE, and the explained 
phenotypic variance of each QTL per environment. In addition 
to determining phenotypic plasticity as QxE, a second method to 
determine plasticity QTL was used as described by (Tetard- Jones 
et al. , 201 1), by QTL mapping the difference in the mean phenotypic 
values per line between treatments. 



Confirming reproducibility of four QTL clusters 

To confirm the reproducibility of the major QTL detected in the 
full RIL screening experiment, this work selected four QTL clus- 
ters: on chromosome 3 at 38^42 cM, on chromosome 7 at 30^40 
cM, on chromosome 8 at 85-95 cM, and on chromosome 9 at 70-84 
cM. The whole population was genetically classified into 16 groups 
based on all possible allelic combinations at the four selected QTL. 
Thereafter, for every tested QTL, phenotypic data of RILs with 
contrasting genotypes for one QTL, but similar genotypes for the 
other QTL, were compared. For example, to test for the QTL on 
chromosome 3, ANNN RILs were compared with BNNN RILs in 
paired groups, so AAAA with BAAA, ABAA with BBAA, ABBA 
with BBBA, etc. The 27 RILs with the highest and lowest average 
values at each tested QTL were selected and grown as described. 
For all measured traits, a correlation analysis between traits meas- 
ured in the control environments and traits measured in the drought 
environments of the full RIL screening experiment and the QTL 
reproducibility experiments was used to test for a significantly 
similar response to the treatment as a confirmation of the level of 
reproducibility. 

Results 

Phenotyping the RIL population 

The results obtained from the pilot experiment (data not 
shown) indicated there was ample phenotypic variation for 
drought response, which justified phenotyping the whole RIL 
population. A total of 20 traits related to growth and per- 
formance of plants were analysed under control and drought 
environments. Fig. 1 shows the frequency distributions of 



the measured traits over the whole population. Transgression 
beyond both parental lines was observed for most of the traits 
except for root volume, RDW, HDW, SDW, S/SR, and LWC, 
where transgression was only in one direction. 

The drought treatment decreased fresh weight, leaf num- 
ber, leaf area, LDW, root length, and stomatal conductance 
and increased R/S (Fig. 1, Table 1). For stomatal conduct- 
ance, the reduction in the L58 parent was minor and not sig- 
nificant, as was also the case for some of the RILs. 

Correlation analysis of all measured traits in this experi- 
ment was performed to unveil the genetic and physiological 
relationships of the various traits (Table 2). The correla- 
tions may exist because of similar physiological mechanisms 
or pleiotropy; however, correlations can also be caused by 
genetic linkage of loci affecting different traits, which are 
not physiologically related or pleiotropic. For instance, the 
analysis showed that, in the control environment, chlorophyll 
content was positively correlated with root diameter, which is 
hard to envision being because of pleiotropy. 

The correlation observed between root length and S/SR 
was positive in the control environment (longer roots contrib- 
uting to relatively more shoots) but negative in the drought 
environment, indicating a proportionally higher investment 
in roots. Under both environments, LWC was negatively cor- 
related with SDW, root length, root volume and RDW, while 
it was positively correlated with LFW and negatively corre- 
lated with stomatal conductance under drought conditions. 
Stomatal conductance was negatively correlated with LDW 
in the drought environment, but positively correlated under 
control conditions. In general, plants with longer root sys- 
tems had higher plant DW. 

As expected, all traits measured in control and drought 
environments showed a positive correlation, except for LWC. 

Mapping QTL with main effects and QxE 

In total 54 QTL were mapped for the traits analysed under 
control and drought environments (Table 3, Fig. 2). Six QTL — 
STC1, LA4, SLA1, RD3, RLISDW3, and SISR2— had oppo- 
site allelic effects when comparing both environments. The 
phenotypic effects of three QTL — LAI, SLA1, and S/SR2 — 
were 9-, 101-, and 15-times higher, respectively, in one environ- 
ment than the other (Table 3). SLA1 colocated with 3.4DW1, 
with the alleles increasing the trait values in the control envi- 
ronment from L58 and R-o-18, respectively. Four QTL were 
mapped for chlorophyll content, of which CHL1, CHL2, and 
CHL3 showed the highest effect from the L58 allele, while for 
CHL4 the R-o-18 allele had the highest effect in both environ- 
ments. Hypocotyl length was mapped to four loci, with the 
R-o-18 alleles contributing most to increased hypocotyl length. 
In total, 1 1 QTL clusters were observed, of which seven com- 
prised at least three colocating QTL (Table 3, Fig. 2). 

Stomatal conductance QTL (STC1 ) and fitness trade 
offs in the drought environment 

The correlation analysis showed that, in the control environ- 
ment, stomatal conductance was positively correlated with 
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Fig. 1. Frequency distributions of the non-normalized trait values for the L58 x R-o-1 8 recombinant inbred line population under control 
(C) (dark grey), and drought (D) (light grey) conditions. Vertical axes indicate the number of lines per trait value class, and horizontal axes 
indicate the different trait value classes. L, L58; R, R-o-1 8. Trait abbreviations are given with Table 1 . 



LDW and SDW. On the other hand, in the drought environ- associated with altering the trait-value-enhancing allele for 
ment, stomatal conductance showed a significant negative STC1 from R-o-1 8 in the control environment to L58 in the 
correlation with LDW and SDW. These correlations were drought environment. The trait-value -enhancing alleles for 
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Table 1 . Parental line means for the analysed traits and performance of the phenotyped RIL population under control and drought 
conditions 

For all traits, three replicates were measured. Min and Max indicate the lowest and highest RIL values; mean and standard deviation (SD) is for 
all RILs; h 2 indicates broad sense heritability. S, significant; NS, nonsignificant; -, not applicable (ANOVA and heritability cannot be calculated 
because averages and not replications were used). CHL, chlorophyll content; HDW, hypocotyl dry weight; HL, hypocotyl length; LA, leaf 
area; LDW, leaf dry weight; LFW, leaf fresh weight; LN, leaf number; LWC, leaf water content; RD, root diameter; RDW, root dry weight; RIL, 
recombinant inbred line; RL/SDW, root length-to-shoot dry weight ratio; RL, root system length; RV, root volume; R/S, root-to-shoot dry weight 
ratio; SDW, shoot dry weight; S/SR, shoot-to-plant dry weight ratio; SLA, specific leaf area; STC, stomatal conductance; DWR, 3,4DW vs. LDW 
ratio; 3,4DW, 3rd and 4th leaf dry weight. 
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the QTL colocating with STC1 (LFW2, LDW2, 3,4DW2) 
were R-o-18 in both environments (Table 3, Fig. 2). This 
means that the R-o-18 alleles for these loci were enhancing 
fitness both under control and drought conditions, although 
having contrasting phenotypic effects on stomatal conduct- 
ance when comparing both conditions. 



Mapping QTL underlying plasticity 

Seventeen of the mapped QTL showed a significant QxE 
effect (Table 3) indicating the loci contributing to pheno- 
typic plasticity between both environments. In addition to 
the GenStat method to determine these plasticity loci, an 
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Table 3. QTL detected in the L58 x R-o-18 RIL population for the traits described in Table 1, using the multienvironment analysis approach 

Per trait, QTL are numbered according to chromosome. R 2 is the percentage of total phenotypic variance explained by each QTL. Effects 
with positive values represent a positive contribution of the R-o-1 8 allele to the trait value and those with negative values represent a positive 
contribution of the L58 allele to the trait value. Highlighted results show significant QxE effects. Ratio refers to the ratios between the effects of 
each QTL in both environments. Trait abbreviations are given with Table 1 . 



Trait 



QTL 



Control 



Drought 



Ratio 



Name 



Linkage 
group 



Position of 
highest peak (cM) 



-loglOP 



Effect 



Effect 



CHL 



STC 
LN 

LA 



LFW 



LDW 



3.4DW 



DWR 



SDW 



SLA 
HL 



HDW 

RL 

RV 
RD 

RL/SDW 

RDW 

R/S 

S/SR 



CHL1 
CHL2 
CHL3 
ZHL4 




-0.855 
-0.885 
-1 .304 
1.011 
-9.754 
11.103 
4.650 
-0.358 
-1.251 
1.321 
-0.394 
0.062 
0.092 
0.120 
0.010 
0.016 
0.023 
0.003 
-0.01 1 
0.007 
0.024 
0.020 
-0.076 
0.011 
0.016 
0.025 
0.201 
0.245 
0.176 
0.265 
0.200 
0.003 
0.002 
0.002 
0.003 
30.215 
32.583 
62.124 
0.056 
0.124 
-0.010 
-0.008 
0.010 
-0.007 
-133.0 
-236.8 
-142.8 
182.5 
0.005 
0.006 
-0.014 
-0.021 
-0.010 
-0.013 



5.4 
5.6 

13.4 
9.7 
7.9 

23.5 
8.1 
0.9 

10.7 

11. 
1.1 
3.7 
8.1 

13.7 
3.0 
8.4 

16.5 
1.3 

14.9 
5.5 
4.2 
3.0 

42.8 
2.9 
7.0 

16.3 
0.0 

12.5 
6.5 

14.7 
8.3 



6.8 
5.9 

11.3 
2.7 
3.2 

11.5 
2.4 

11 
4.6 
3.3 
5.0 
2.6 
3.1 

10. 
3.6 
5.9 
6.7 

10.3 
3.7 
8.5 
4.7 
8.2 



5.0 



3.7 



3.6 



-101.4 




-14.S 
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Fig. 2. A clustered heat map showing the -log10P profiles of the measured traits. Columns indicate the 10 chromosomes in 
centiMorgans, ascending from left to right; rows indicate individual trait -log10P profiles. A colour scale is used to indicate the QTL 
significance corresponding to the -log10P score: red and black represent a positive effect on the trait value from the R-o-18 allele; 
blue and green represent a positive effect on the trait value from the L58 allele. Bar width indicates the significance interval of the QTL. 
Hierarchical clustering, shown on the left, reflects the correlation between traits based on the QTL profiles. Right dotted line of heavier 
weight in A7 indicates the QTL with an antagonistic fitness effect; the other five dotted lines refer to QTLs confirmed in the reproducibility 
experiment. C and D refer to control and drought environments, respectively. Trait abbreviations are given with Table 1 . 



alternative method to describe QTL that are affected by the 
environments was suggested by Tetard-Jones et al. (2011). 
This uses the differences between the trait-value averages of 
the lines in the two environments to determine QTL. Using 
this procedure, 15 plasticity QTL were mapped (Table 4), with 
nine of them colocating with previously mapped QTL, six of 
which were found to show QxE (Table 3). Thus, this analysis 
detected six new plasticity QTL, which did not exceed the sta- 
tistical significance levels with the GenStat method. 

Reproducibility 

From the 11 QTL clusters that were mapped, four were 
selected to be tested for reproducibility in a subsequent exper- 
iment. The first cluster mapped to A3, including the RD1 and 
RLISDW1 QTL for both positively correlated traits, with 
trait-value-enhancing effects from the L58 alleles. Moreover, 
LFW1, LDW1, and SDW1, all contributing to shoot bio- 
mass, were mapped to the same cluster, with positive alleles 



coming from R-o-18. The second cluster was mapped to 
A7— composed of LN1, LA2, 3.4DW2, DWR3, RJS1, and 
RLISDW2 — all with a positive contribution of the R-o-18 
allele except for LN1. This is in line with the negative correla- 
tion of leaf number with the other traits. The third cluster, on 
A8, included eight colocating QTL — LFW3, LDW3, SDW3, 
LA3, RL3, RV12, RD4,and RDW2—oi which the RD4 L58 
allele increased the trait value, while for the other QTL, the 
R-o-18 allele increased the trait value, in line with the nega- 
tive correlation of RD with the other traits. The fourth clus- 
ter included three QTL — CHL3, RJS2, and S/SR2— mapping 
to A9. The S/SR ratio showed a negative correlation between 
control and drought environments and therefore the trait- 
value-enhancing effect of S/SR2 in the drought environment 
came from the L58 allele, whereas in the control environment 
it came from the R-o-18 allele. 

In total, 27 lines were selected from the RIL population 
to properly represent the 16 possible genotypes for all allelic 
combinations for the four selected QTL clusters. These lines 
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Table 4. QTL mapped for phenotypic plasticity in the L58 x R-o-18 RIL population 

Plasticity was calculated as described by Tetard-Jones et al. (201 1), as the difference in the mean phenotype between different treatments per 
trait. QTLs are numbered according to chromosome. R 2 is the percentage of total plastic variance explained by each QTL. Effects with positive 
values represent a positive contribution of the R-o-1 8 allele to the trait value and those with negative values represent a positive contribution of 
the L58 allele to the trait value. Highlighted QTL were mapped before using the multienvironment analysis approach (Table 3). Trait abbreviations 
are given with Table 1 . 



Trait 



QTL name 



Locus 



Chromosome 



Position (cM) 



Effect 



-loglOP 



STC 

LA 

LDW 

3.4DW 

SDW 

SLA 

HL 
RL 



RV 
RD 

RL/SDW 
R/S 




903607|9917837 
E3835M1 1 
E3850M9 
Ra2A01 -A7 
E3850M9 
E3749M6 
BrlD101239-A5 
902225|9924661 
E3732M5 
E3850M9 
E3416M22 
E3732M5 
89901 5|991 8455 
C7P119 
BrlD10177-A9 




Table 5. Correlation analysis between the control conditions and 
the drought conditions of the full RIL screening experiment and 
the reproducibility experiment 

Highlighted results refer to significant correlations: dark grey, P < 0.01 ; 
light grey, P < 0.05. Trait abbreviations are given with Table 1 . 



Trait 


Control 


Drought 


STC 




0.494 


LN 


0.852 


0.758 


LA 




0.220 


LFW 


0.122 


-0.027 


LDW 




0.258 


3.4DW 


0.474 


0.239 


DWR 


0.296 


0.128 


SDW 




0.259 


SLA 


0.496 


0.159 


HL 


0.850 


0.847 


HDW 


0.211 


0.324 


RL 


0.209 


0.085 


RV 


-0.002 


0.292 


RD 


0.071 


-0.107 


RL7SDW 


0.095 


0.378 


RDW 




0.251 


R/S 


0.310 


S/SR 




LWC 


0.105 


0.088 



were regrown under similar conditions and rephenotyped 
(Supplementary Table 1, available at JXB online). A corre- 
lation analysis (Table 5) between traits measured in the two 
control environments and between traits measured in the two 
drought environments of the full RIL screening and QTL 



reproducibility experiments showed that all traits were posi- 
tively correlated in at least one environment, but often both, 
except for fresh weight, LWC, and root diameter. This indi- 
cates that the phenotyping was robust and the detected QTL 
clusters are reproducible, making them attractive candidates 
for further gene cloning experiments. 

Discussion 

The current study was carried out in a greenhouse using 
pots filled with sand. This type of pot experiment is a rea- 
sonable compromise to avoid the difficulty of phenotyping 
roots in natural field environments and the unnatural con- 
ditions present in hydroponics, aeroponics, or agar plates 
(Tuberosa, 2012). However, aspects of root growth in this pot 
system would have still been substantially different from field 
conditions. 

Upon screening the RIL population, this work found sig- 
nificant GxE between control and drought environments for 
stomatal conductance, leaf number, root length, and root 
diameter. This GxE was reflected in QxE detected using 
the MEA approach for these traits, except for leaf number. 
MEA is more powerful than the traditional single environ- 
ment analysis in detecting more significant QTL with higher 
explained variance. An additional advantage is that it allows 
quantification of QxE, because it accounts for GxE and 
tests all detected QTL in all environments and thus shows 
their effects in each environment (Crossa and Federer, 2012). 
QxE occurs if the QTL effects are strongly expressed in one 
environment but weakly in another, or if the QTL has oppo- 
site effects on the same trait in two different environments 
(Mackay, 2001; MacMillan et al, 2006; Zhang et al, 2010). 
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Examples of the first case are LFW1, LAI, LDW1, 3.4DW1, 
RL1, RL2, RV1, RIS2, HDW3, and SDW1, while examples 
of the latter case are found for LA4, SLA1, RD3, and S/SR2. 
The latter kind of QxE obstructs the transferability of QTL 
mapping results from one environment to another (Mackay, 
2001), as selection will be in opposite directions in the two 
environments. 

Knowing about the QTL with opposite effects on several 
traits in different environments, also known as antagonistic 
pleiotropy, is of great importance in breeding programmes 
because breeding for one trait might negatively affect other 
traits (Rose, 1982; Juenger, 2013). The QTL cluster mapped 
at the bottom of A7 included a stomatal conductance QTL 
(STC1), which showed signs of antagonistic pleiotropy, with 
the R-o-18 allele increasing stomatal conductance under con- 
trol conditions and decreasing it under drought conditions, 
while having positive effects on biomass under both environ- 
ments through the colocated LDW3, SDW2, and 3.4DW3 
QTL. However, the similar effect on biomass and the con- 
trasting effect on stomatal conductance could also mean 
these traits are not allelic, but the result of close linkage of 
two loci. Further analysis should reveal this. 

Stomatal conductance showed clear plasticity, decreas- 
ing significantly in the drought environment. Such response 
is generally correlated with reduced photosynthesis but also 
with reduced water loss as an adaptive response to drought 
(Chaves et al, 2003; Condon et al, 2004; Tardieu, 2013). Due 
to the colocation or antagonistic pleiotropy of the shoot bio- 
mass QTL with STC1, when comparing both environments, 
stomatal conductance was negatively correlated with shoot 
biomass (Table 2B). This reflects an interesting fitness advan- 
tage for plants carrying the R-o-18 allele at this QTL cluster, 
meaning that under drought conditions they show relatively 
reduced stomatal conductance (contributing to increased 
drought tolerance) accompanied with relatively increased 
shoot biomass, compared to plants carrying the L58 allele. 

Recently, the plasticity and the evolution of flowering 
time and water use efficiency (WUE) has been investigated 
in B. rapa under drought environments (Franks, 201 1), and 
the relationship between circadian rhythm, vegetative, and 
reproductive traits, and leaf gas exchange with the varia- 
tion of WUE in different watering regimes has been inves- 
tigated (Edwards et al, 2012). The negative correlation 
that this work found for stomatal conductance and shoot 
biomass under drought was also observed by (Edwards 
et al, 2012), although this was not significant in their 
study. It also agrees with the positive correlation between 
WUE and biomass in the drought environment found by 
these authors and the colocation of WUE and stomatal 
conductance QTL mapped in B. rapa grown under warm 
and long-day conditions (Edwards et al, 2011). Although 
the preferred targets for crop improvement in marker- 
assisted breeding are generally constitutively expressed 
QTL (Bernardo, 2008), this QTL cluster is attractive to 
select for, even if it is not constitutive in view of the QxE 
observed for STC1, because the allele from R-o-18 con- 
tributes to increased drought tolerance without having fit- 
ness costs due to reducing biomass. 



The leaf area response and the underlying QTL in both 
environments were confirmed by the positive correlation 
observed between the full RIL screening and QTL reproduc- 
ibility experiments. Stomatal closure and limited expansion 
of young leaves under drought have an indirect negative effect 
on root growth (Chaves et al, 2003; Roycewicz and Malamy, 
2012). This was observed by the reduction in root length, 
concomitant with an increase in root diameter in the drought 
treatment, corresponding to similar observations reported 
before for Brassica and other crops (Zhu et al. , 201 1 ; Edwards 
et al, 2012; Poorter et al, 2012). It thus appears that, under 
drought stress in pots, B. rapa does not invest in longer roots 
to take up more water, but in thicker roots to act as a water 
storage buffer. 

Under drought, the R/S ratio increased compared to the 
well-watered conditions. Biomass allocation under limiting 
environments can be explained by a functional biomass equi- 
librium when plants allocate more biomass to roots when the 
factor limiting growth is below ground (e.g. water or nutri- 
ent shortage), to enhance the uptake of that limiting factor 
(Poorter et al, 2012). The correlation of the R/S ratio with 
drought tolerance has previously also been documented for 
Arabidopsis and tobacco (Werner et al, 2010), as well as 
B. rapa (Kage et al, 2004; Edwards et al, 2012). 

Of the traits examined, GxE was found for most of them, 
either in the full RIL screening or the reproducibility experi- 
ment (Table 1). With so many traits for which GxE was found, 
it is not surprising that QTL with QxE were also found, indi- 
cating plasticity for many traits. This work used two ways to 
detect QTL related to phenotypic plasticity, first using the 
MEA approach (Table 3) and subsequently using the differ- 
ence between average values per line when comparing both 
treatments per lines (Table 4). As previously found by Tetard- 
Jones et al. (201 1), there is considerable overlap between both 
methods, but the latter method also detects some novel QTL 
not found previously. This is probably due to the additional 
statistical power that can be gained by directly using the phe- 
notypic difference values for mapping, meaning that QTL 
that did not exceed the threshold in the MEA approach will 
be detected. 

Although almost all traits showed a positive correlation 
between the results from the full RIL screening experiment 
and the reproducibility experiment, confirming the initial 
results, this was not the case for leaf fresh weight and root 
diameter, suggesting a high level of GxE for those traits, or 
for water content, where a high environmental effect probably 
prevented mapping a QTL for this trait. The only trait for 
which no plasticity QTL was found was chlorophyll content 
(CHL), which was in line with the inability to detect GxE for 
this trait (Table 1). However, there was genetic variation for 
CHL, with four detected QTL (Table 3). There is also a differ- 
ence in CHL between drought and control conditions, which 
agrees with previous observations for four Brassica species 
(Ashraf and Mehmood, 1990), but the genotypes appeared 
to respond similarly to drought exposure by decreasing CHL, 
explaining the lack of GxE. The four CHL QTL mapped 
to regions previously identified to contain QTL for chloro- 
phyll a and b content in B. rapa (Ge et al, 2012), with CHL1 
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colocating with one of the three QTL previously identified for 
chlorophyll fluorescence (Edwards et al, 201 1). 

Increasing crop productivity under drought conditions is 
the ultimate goal for marker-assisted breeding programmes. 
In that respect, the significant antagonistic effect of relatively 
reduced stomatal conductance along with relatively higher 
shoot biomass under drought conditions due to the STC1I 
shoot biomass locus at the bottom of chromosome A7 is 
very interesting, as it suggests that selection on reduced water 
loss during drought, through reduced stomatal transpiration, 
is expected to have disproportionately little effect on shoot 
biomass reduction, which is a favourable combination. In 
addition, this work reported many QTL underlying several 
morphological and physiological traits, which appeared to 
be robust and thus provided the first step towards identify- 
ing genes governing those traits. The availability of the whole 
B. rapa genome sequence (Wang et al, 2011) together with 
possible comparative alignment with the related model spe- 
cies A. thaliana (Schranz et al. , 2006) will facilitate fine map- 
ping and cloning of candidate genes underlying the desired 
QTL. This approach will not only be useful in breeding 
B, rapa, but also in breeding other closely related species like 
B. juncea and B. napus (Cheng et al, 201 1; Li et al, 2013). 

Supplementary material 

Supplementary data are available at JXB online. 

Supplementary Table SI. Phenotypic data of the selected 
lines for the analysed traits by full RIL screening and QTL 
reproducibility under control and drought conditions. 
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